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Non linear excess conductivity of Bi 2 Sr 2 Ca n _iCu n 02 n +4+x (n = 1, 2) thin films 
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The suppression of excess conductivity with electric field is studied for Bi2Sr2Ca„_iCu n 02 n +4+a; 
(n = 1, 2) thin films. A pulse-probe technique is used, which allows for an estimate of the sample 
temperature. The characteristic electric field for fluctuations suppression is found well below the 
expected value for all samples. For the n = 1 material, a scaling of the excess conductivity with 
electric field and temperature is obtained, similar to the scaling under strong magnetic field. 
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I. INTRODUCTION 

High-T c superconductors exhibit strong superconduct- 
ing fluctuations at the superconducting temperature, as 
a result of the high ratio of the transition temperature to 
the condensation energy - essentially due to their short 
coherence length. Such fluctuations have been exten- 
sively studied in the weak and strong fluctuation lim- 
its (respectively the Gaussian and the critical regime). 
There is now a rising interest in the possibility to sup- 
press superconductivity by the application of an electric 
field, and the study of the electric field dependence of 
the fluctuations appears as a first step towards this goal. 
Early theoretical work has focused on the reduction of 
the excess conductivity in the Gaussian regime, by the 
application of a transport current approaching the de- 
pairing current (non linear conductivity). The electric 
field dependence of the Aslamazov-Larkin termpj was 
obtained in Refs. @,0iE1O1> an d the one of the additional 
'anomalous' Maki-Thompson termjg in Ref. 0. Experi- 
mental results on Al films demonstrated the suppression 
of fluctuations by the electric field in the isotropic two- 
dimensional limit 0,0. The theory, based on a Langevin 
equation for the order parameter, was extended to below 
T c in Ref. |lfj|. Later, several works considered the con- 
tribution of the interaction of fluctuations in the critical 
region and, recently, the expression for the fluctuation 
conductivity in the transition region was revisited and 
the similarity of the results with those obtained for the 
broadening of the transition under a strong magnetic field 
was pointed out (see Ref. an d refs therein). Such 
experiments on high-T c superconductors are difficult, as 
they require high current densities, which generally in- 
duce an uncontrolled heating of the sample. As a result, 
only few experimental works have investigated the effect 
of the electrical field on the resistive transition of the 
cuprates. fn Ref. 01 > a field- induced crossover between 
a three-dimensional and two-dimensional behavior was 
reported for a YBa2Cu307_a; single crystal, in the tem- 
perature range e = (T — T c )/T c < 2 1CT 2 . fn Ref. P3, 
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the non-linear resistivity in the range 0<e<710 -2 
was used to obtained the characteristic depairing field 
of L^Si^CaC^Og+a; single crystal, ft was found to be 
about fO times smaller than the expected value and its 
temperature dependence was clearly different from the 
e 3 / 2 dependence predicted in Ref. 

In this contribution, we compare the excess con- 
ductivity suppression for this system with the one in 
Bi2Sr2CuOg+a:, where the investigation of temperature 
intervals as large as —0.2 < e < 0.6 was possible due to 
an electrical critical field more than one order of magni- 
tude smaller. 



II. EXPERIMENT 

The measurements were carried out on 
Bi 2 Sr 2 Ca„_iCu„02„+4+x (n = 1, 2) epitaxial, c- 
axis oriented thin films. They were grown on heated 
SrTi03 substrates, by reactive rf sputtering with an 
oxygen rich plasma (Ref. and refs therein). After de- 
position of typically 2500 A thick films, Au contacts were 
sputtered, and the samples were patterned in the four 
contact transport geometry, with a current carrying strip 
of typical width and length (L) of 80 /zm and 200 /j,m. 
The annealing of the samples at low temperature was 
used to set the doping state, which was determined from 
the normal state resistivity temperature dependence|l5|. 
The thin film with n = 2 (sample II) was overdoped. 
A thin film with n — 1 was studied for two different 
doping states: an overdoped one (sample la) and an 
underdoped one (sample lb). Two other samples with 
n — 1 (samples Ic and Id), respectively close to optimal 
doping and strongly overdoped, were also investigated. 
The sample parameters are summarized in Table [i] 

The resistivity in the limit of vanishing current den- 
sity was obtained using a 10 /iA ac current and a lock-in 
amplifier detection. The investigation for large current 
densities was carried out using the pulsed current tech- 
nique. The current pulses where 10 fis long, with a repe- 
tition rate r = 10 and the maximum current value was 
100 mA. Both the voltage on a reference resistor fed with 
the measuring current and the one on the sample, after 
amplification by home-made amplifiers, were recorded, 
using a 2.5 MHz, 16 bits digital acquisition card. The 
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measurement of the resistance, obtained from the ratio 
of the two voltages, was linear within 0.3% over the cur- 
rent range. The non linearity was corrected from the raw 
measurements. Evaluating the temperature increase due 
to the measuring current is crucial for the high current 
density transport experiments, as an apparent decrease 
of the excess conductivity with current may occur, due to 
the sample temperature increase only. A reliable deter- 
mination of the non linear excess conductivity requires a 
method for measuring the sample temperature. For this 
purpose, we modified the technique to a 'pulse-probe' 
one: after the completion of the main 10 /is pulse with 
current density J and the measurement at the end of this 
pulse, a probe pulse with current J s <C J is performed 
immediately after. A typical secondary current pulse was 
J s ~ J/50. Neglecting the sample heating by current J s , 
the value of the resistance during the secondary pulse 
provides a measure for the sample temperature about 
1 /is after the main pulse, using the ac resistivity data 
as a thermometer. Such a method does not, however, 
catch the fast temperature relaxation which occurs be- 
tween the main and the probe pulses. This relaxation, 
with a typical relaxation time in the ns range. 16], may 
be estimated as AT ~ Ri 2 /AX, where A is the sample 
surface and A ~ 10 7 WK _1 m -2 is the boundary conduc- 
tance, yielding AT ~ 0.01 K and 0.1 K for samples with 
n = 1 and n = 2 respectively. 



III. RESULTS AND DISCUSSION 

The temperature increase for the over-doped 
Bi2Sr2CaCu208+ K thin film (sample II) is shown 
in Fig. \%)p. The resistance of the current contacts 
for this sample was 20 il. As can be seen from the 
comparison with the resistance of the sample and, also, 
from the scaling of the sample heating with J 2 , the 
data may be described using the simple evaluation: 
(^sample ^experiment) C< R{T Q ,J) J 2 below - 100 
K. This demonstrates that the strip temperature is 
essentially determined by the heating power in the 
film constriction rather than by the one dissipated in 
the current contacts. Using the probe value as the 
actual sample temperature, a slight broadening of the 
transition remains which must be attributed to the 
transport current (Fig. |2Jd). 

The reduced excess conductivity was obtained as : 
a'(J)/a'(J = 0) = R(Q)(Rn - R(J))/R(J)(R n - R(0j), 
where R n is the normal state resistance. R n was ob- 
tained from the extrapolation to below T = 130 K of 
the empirical function R = a + b T a for this overdoped 
sample 15], which fitted the data in the temperature in- 
terval 130-300 K with a = 1.13. As can be seen in Fig. [31 
the broadening of the transition corresponds to a drop 
of the reduced excess conductivity from the normal state 
asymptotics a' (J) / a' (0) — 1 to zero with decreasing tem- 
perature. A similar, but far more dramatic, behavior is 
observed for the F^S^CuOg+i thin films. Samples la 



and Id were overdoped, with a = 1.08 and a = 1.10 
respectively. Samples lb and Ic were underdoped. For 
these samples, the normal state resistivity was fitted us- 
ing an activated law. Indeed, as in the case of under- 
doped YBa 2 Cu 3 6+a; [l3 and Bi 2 Sr 2 CaCu20 8+a; [l3, it 
was found that the resistivity is well described by the 
activated law po + b exp(— A/T) below T*, and by a 
linear temperature dependence above. T* ~ 65 K and 
T* ~ 85 K was found respectively for samples lb and Ic. 
As in Ref. [13 , it was found that A ~ T* . For all sam- 
ples with n — 1, the sample heating was found to increase 
with decreasing temperature (Fig-QJi): this is likely due 
to the decrease of the sample specific heat and of the ther- 
mal conductivity at low temperature. For these samples, 
the resistivity at temperatures well above the mid-point 
transition is clearly affected by the transport current, 
which points to the suppression of the excess conductiv- 
ity by the current. Finally, we stress that the resistivity 
curves obtained at constant J are not equivalent to the 
constant E curves usually obtained from theoretical com- 
putations, and that the electric field [E = R I / L) must 
be calculated for each data point in order to compare 
with the theoretical results. 

Before we compare these results with theory, let us ex- 
amine the experimental artifacts which may affect the 
results given above. First, the normal state resistiv- 
ity, as extrapolated from the high temperature values 
to the superconducting transition, may not be exact. It 
can be shown that the error made in the reduced excess 
conductivity is of the order of (8R n / a') 2 . So, the er- 
ror which results from the normal state indetermination 
drops strongly as one enters the superconducting tran- 
sition, as evidenced by the error bars in Fig. [31 Also, 
despite very different normal state approximations were 
used (power law or activated ones), the characteristic 
field obtained from the excess conductivity was invari- 
ably well below the theoretical one. Thus, it is un- 
likely that the uncertainty on the normal state resistivity 
value can account for this central result of our measure- 
ments. Then, the accuracy of the sample temperature 
may be questioned. The error made in the computation 
of the temperature after the main pulse is critical only 
for sample II. At the middle of the transition for this 
sample, the broadening at the highest current density is 
about 0.6 K. This is larger than the temperature increase 
(0.2 K)obtained from the probe signal and, due to the 
high sensitivity of the method in the transition region, 
noticeably larger than the expected error on the sam- 
ple temperature. Then, the fast temperature relaxation 
that occurs after the main pulse should be estimated. 
After correcting the temperature, the sample resistance 
values above the superconducting transition (say, 14 K 
and 100 K for sample la and II respectively ) where the 
effect of current reduces to that of heating, agree within 
a temperature shift of about 0.2 K. While such a tem- 
perature difference would yield a sizeable correction of 
the results for sample II, this is clearly not so for sample 
with n — 1. Finally, it may be argued that the sample 
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inhomogeneities affect the transition and its broadening 
with electric field. The finite width of the transition ( 
~ 8 K, ~ 2 K, ~ 3 K, ~ 3 K and ~ 5 K for the 10%- 
90% transition completion of samples II, la, Id, lb and Ic 
respectively) could be partly attributed to the inhomo- 
geneities that are present in the sample. While it should 
be stressed that, a contrario, a narrow transition does 
not necessarily mean a more homogeneous sample and 
that intrinsic fluctuations also result in a broadening of 
the transition (see Fig. Eh), this could be indeed a serious 
limitation in the present case. As will be seen later, the 
universal behavior observed for the n = 1 samples, inde- 
pendently of their oxygen concentration, rules out such 
an interpretation in this case. 

The theoretical results considering only non interacting 
fluctuations, i.e. Gaussian ones, for a 2D system were 
given in Refs. 0,0, 0,13- I n this case, the normalized 
excess conductivity is simply given by : 

/>oo 

a'(T,E)/a'{T,0) = / dxexp{-x - [E / E c (T)] 2 x 3 } 
Jo 

(!) 

where E C (T) = E e 3 ' 2 , E = 16V3k B T c /ir e £ 
and e = (T — T c )/T c . The Gaussian result is found to 
account roughly for the data for all samples (Fig. |3J). 
However, the characteristic field obtained from these 
fits is Eq ~ 10 3 Vm" 1 for samples with n = 1 and 
Eq ~ 3 10 5 Vm _1 for sample II. These values are con- 
siderably smaller than the theoretical estimates Eq = 
10 6 Vm- 1 and E = 3 10 7 Vm" 1 (using £ = 40 A and 
£o = 20 Arespectively) . As shown in Ref. there is 
an enhanced suppression of the excess conductivity, be- 
low the zero-field transition temperature, when critical 
fluctuations are taken into account within the Hartree 
approximation. Such an enhancement may result in a de- 
crease of the apparent value for Eq when only Gaussian 
fluctuations are considered. Within the Hartree approxi- 
mation, the magnitude of the fluctuation renormalization 
is essentially determined by the Ginzburg number and 
the anisotropy of the superconductivity. These parame- 
ters determine the shift of the superconducting transition 
temperature with respect to the bare mean field temper- 
ature and, associated to this shift, there is an enhance- 
ment of the superconducting fluctuations. For sample 
II, using s = 15.35 A, £ 0c = 0.2 A, k = 100 for the su- 
perconducting plane separation, the transverse coherence 
length and the Ginzburg Landau parameter respectively, 
it is found — as a result of the quasi-2D character for 
this material — that the bare mean field characteristic 
temperature calculated according Ref. 0, 0] is larger 
by about 30 K than the superconducting transition tem- 
perature (in Ref. the fluctuating modes are cut for 
k 2 > c £q an d w e have used here c = 1). Such a differ- 
ence results in a large enhancement of the superconduct- 
ing fluctuations with respect to the Gaussian approxima- 
tion. This clearly fails to account for the experimental 
data taken at E — (Fig. [2b)- Critical fluctuations from 
the Hartree approximation and the raw superconducting 
parameters given above are then clearly strongly over- 



estimated in the region of interest for our data, and it 
is not surprising that we were unable to obtain a satis- 
fying fit of the reduced excess conductivity in this case. 
However, the Hartree approximation results were used 
in Ref. [19) to provide convincing fits of the resistivity 
data under a magnetic field. The anisotropy parameter 
obtained from these fits was r = (2£ i q c /s) 2 = 5 10 -3 and 
the Ginzburg number (using the notation from Ref. |ll| ) 
was g = G 2 /AT C = 4.5 10~ 5 K -1 , while, from the raw 
superconducting parameters given above, one obtains 
r = 1.7 10~ 4 and g = 6.7 10~ 4 K _1 . These differ- 
ences both contribute to a decrease of the magnitude 
of the superconducting fluctuations, and the prediction 
for the zero field resistivity using the same parameter 
values as Ref. [l9j now agrees more reasonably with 
the experimental data (Fig. |2Jd) (although some part of 
the transition broadening is likely due to the sample 
inhomogeneity) . Following the formalism of Ref. [llj . 
the reduced excess conductivity for such r and g val- 
ues is now essentially that of the Gaussian regime (here, 
E < 5 10 3 Vm _1 ) and the experimental data may be 
fitted using E Q (T C ) = 4.7 10 5 Vm" 1 (Fig. Eft). 

In the case of samples with n = 1, using the follow- 
ing parameters : s = 12.3 A, £o c = 2 A, k — 100, the 
bare mean field transition temperature is only shifted by 
about 0.4 K with respect to the superconducting tran- 
sition temperature. The critical region should then be 
much narrower than for sample II and, logically, the 
Gaussian and the Hartree approximations give close re- 
sults for Eq (Fig. |3Ji), i.e. well below the above theoreti- 
cal estimate for all n — 1 samples (Table QJ. The Results 
of Ref. [ll| should also be valid also for T < T c where 
the normalized excess conductivity is zero. As pointed 
out in Ref. [ll| , the treatment of the electric field depen- 
dence of the excess conductivity within the Hartree ap- 
proximation naturally yields results comparable to those 
obtained in the case of a magnetic field, within the same 
approximation^^. As a result, one may expect that 
the excess conductivity also scales with temperature and 
electric field, as was observed[20j in the case of a mag- 
netic field. As shown in Fig. the fluctuation conduc- 
tivit y is found to scale according to the two-dimensional 
law [HI: a' oc {T/E) 1 / 2 ^ - T C (E))/(T E) 1 / 2 }, pro- 
vided one takes dT c (E)/dE a constant as given in Ta- 
bled A three-dimensional law (a' oc (T 2 /S) 1/3 ^ r [(T - 
T C (E))/(T E) 2 / 3 }), yields a slightly less satisfying scaling, 
but comparable dT c (E) / dE value. The critical tempera- 
ture field dependence obtained in this way confirms the 
fit of the reduced excess conductivity with the results in 
ref.HU, as the critical field T c0 (dT c / 'dE)^ 1 is also well 
below the theoretical estimate (Table HJ. 

Thus, there are converging experimental pieces of evi- 
dence that the apparent characteristic field — at least for 
samples with n = 1, for which the transition broadening 
is much larger than the error possibly made for the sam- 
ple temperature — is well below the one expected from 
a simple estimate of the depairing field. Several expla- 
nations may be considered for this. First, it is known 
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that the existence of 'hot spots', where the sample tem- 
perature locally exceeds T c , may drive the sample into 
the normal state. Such instabilities must be ruled out 
above T c , i.e. in the analysis of the reduced conduc- 
tivity. Below T c , the characteristic current density for 
which hot spots develop may be estimated to be [2l| 
10 10 (T c — T) 1 / 2 Am~ 2 , well above current densities which 
significantly affect the resistive curve of n = 1 samples. 
Then, as noticed above, the transition width of the sam- 
ple is finite, which may be interpreted as the result of 
a distribution of the superconducting temperature. The 
effect of T c inhomogeneities on the superconducting tran- 
sition is complicated and the full percolation problem 
(with a non-linear conductivity) should be considered in 
the general case. The situation is more simple when the 
resistivity is close to the normal state value, as the as- 
sumption of a uniform current density may be used. Us- 
ing this assumption in the case of the larger current dis- 
played in Fig. and the results of Ref. [Tl| , it is found 
that the resistivity curve for a field E/Eq ~ 10~ 4 , as 
would be expected from the theoretical estimates above, 
is only weakly affected by a distribution of T c with a typi- 
cal width of 1.5 K. In addition, within the hypothesis that 
this T c distribution strongly influences the apparent value 
for the characteristic electric field, one would expect that 
this parameter would be highly sensitive to the oxygen 
doping state. This is clearly not so, as demonstrated in 
the case of the n = 1 samples (Table [|J: although the 
optimal doping might coincide with a maximum in Eq, 
the characteristic electric field is well below the depairing 
estimate for all samples. Another possible explanation 
for the apparent weakness of the critical field is that the 
observed current-induced resistivity is the one of the in- 
tcrgranular material. Such a contribution is unlikely for 
these epitaxial films where grains (typically 0.1 /an large) 
are well oriented by the substrate and exhibit sharp (0°- 
90°) grain boundaries at the atomic level. It is known 
also that granular bulk or thin films materials close to 
T c are invariably in the Ginzburg-Landau regime, where 
the critical current is the one needed for the suppres- 
sion of the order parameter in the grains 22] . Moreover, 
this would require that the intergranular weak links con- 
tribution to the resistivity is comparable to the normal 
state value, which would likely imply a two step transi- 
tion curve (corresponding to the Ginzburg-Landau to the 
Ambegaokar-Baratoff behavior crossover temperature) at 



large currents, which is not observed here. The scaling 
of the data according to the prediction of the depair- 
ing mechanism rather suggests that a microscopic mech- 
anism (at the coherence length scale) should be found. 
Microdomains at this scale or below 23] is one of such 
candidates and it has already been noticed that they may 
affect the magnitude of the depairing current (2^. Along 
this line, it is worth mentioning that the n = 1 material 
also exhibits a peculiar behavior with respect to the su- 
perconducting Nernst effect. Indeed, whereas the Nernst 
effect for YBa2Cu306+ K and Bi 2 Sr 2 CaCu208+ :r shows a 
peak which shifts with magnetic field to lower tempera- 
ture in a way similar to the resistive transition|24|. the 
underdoped n = 1 material has shown a nearly field inde- 
pendent peak[2^|. Moreover, the Nernst effect clearly has 
its maximum well above the resistive transition temper- 
ature (i.e. in the fluctuation regime), whereas this peak 
points — at zero field — towards the completed resistive 
transition in the case of the two former compounds. As 
underlined in Ref. (25|, this discrepancy between the re- 
sistive and the Nernst effect transitions may be due to an 
anomalously large vortex contribution, such as the one 
arising from a vanishing vortex viscosity. Such a mecha- 
nism could account for an enhanced dissipation with cur- 
rent, as observed here. The existence of a phase coher- 
ence length distinct from the amplitude one, as proposed 
in the case of microscopic inhomogeneity 23] , may also 
account for these anomalies. Finally, the contribution of 
the d-wave symmetry to the depairing field should also 
be considered as, in the nodal directions, the depairing 
field is virtually zero: although there is to our knowl- 
edge no theoretical evaluation in this case, we expect a 
decrease of the critical field with respect to the conven- 
tional s-wave symmetry, as well as an in-plane angular 
dependence of the non-linear resistivity. Clearly, further 
studies of the effect of large electric field, in particular 
on compounds which have shown distinct resistive and 
thermodynamic lines — such as Tl2Ba2Cu06+5 26] — are 
needed to clarify these mechanisms. 
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TABLE I: Sample parameters 



Sample 


n 


To(K) 
(K) 


doping 


E a 
Vm' 1 


T c (dT c /dE)- lb 
Vm' 1 


II 


2 


67.4 


overdoped 


4.7 10 b 




la 


1 


5.65 


overdoped 


3.7 10 3 


9.5 10 2 


lb 


1 


9.8 


underdoped c 


7.8 10 3 


1.5 10 3 


Ic 


1 


14.4 


underdoped 


1.5 10 4 


4.8 10 3 


Id 


1 


5.1 


overdoped 


1.8 10 3 


3.8 10 2 




80 100 
T(K) 



FIG. 1: Sample la (a) and sample II (b): temperature differ- 
ence between the sample holder and the film, as measured by 
the probe pulse. 



"from the fit using the results of 111! . 
'from the 2D scaling. 

c same sample as la after vacuum annealing 
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FIG. 2: Sample la (a) and sample II (b): resistivity vs sam- 
ple temperature. Dotted line: the extrapolated normal state 
resistivity a + bT a ; dashed line: Hartree approximation us- 
ing the parameters in Ref. [l9j|: dashed dotted line: Hartree 
approximation using the raw superconducting parameters. 



7 




FIG. 3: Sample la (a) and sample II fb): normalized ex- 
cess conductivity. Dotted line is Ref. 0] Gaussian approxi- 
mation using T c = 5.65 K and Eq = 1900 Vm" 1 for sam- 
ple la and T c = 68 K and E = 3 10 5 Vm" 1 for sam- 
ple II. Full line is ref. [HI Eqs.(24),(29) , using respectively 
E {T C ) = 3700 Vm" 1 and E (T C ) = 4.7 10 5 Vm" 1 . 
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FIG. 4: Sample la: 2D scaling, using linear T C (E). 



